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Abstract 

This paper generalizes some known solitary solutions of a time- 
dependent Hamiltonian in two ways: The time-dependent field can be 
an elliptic function, and the time evolution is obtained for a complete 
set of basis vectors. The latter makes it feasible to consider arbitrary 
initial conditions. The former makes it possible to observe a beating 
caused by the non-linearity of the driving field. 

1 Introduction 

The analytic solutions of Allen and Eberly [1] for the Bioch equations are 
well-known. Similar results for spin-one systems or three-level atoms do 
exist [2] and are derived in terms of the coherence vector of Hioe and Eberly 
[3]- We consider a three- level system with time-dependent external fields 
which enable transitions between two pairs of levels, between (1) and (3), 
respectively between (2) and (3). See the Figure [H This kind of system has 
applications in different domains of physics. Analytic expressions for the time 
evolution of the density matrix are for instance very helpful for understanding 
many of the phenomena observed in light scattering experiments — see for 
instance [3J. In the context of quantum computers the accurate manipulation 
of the state of a quantum system — in this qutrit — is important. 

In the present work the solitary solutions of [2] are generalized in more 
than one way. The external fields are modulated with Jacobi's elliptic func- 
tions. By varying the elliptic modulus k these functions make the bridge 
between periodic functions (cos(ut) and sm(wt)) and single pulses described 
by sech (ut) and/or tanh(W). In addition, a full set of solutions is presented 
instead of just one solution. This makes it possible to take arbitrary initial 
conditions at time t — 0. 
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The next Section presents the time-dependent Hamiltonian and the spe- 
cial solutions. In Section 3, a specific setting is chosen. Section 4 discusses 
the results. The Appendix A contains the explicit expressions which are used 
for the generators of SU(3). The Appendix B explains the method by which 
the special solutions were obtained from known solutions (see the appendix 
of [S]) of the non-linear von Neumann equation 

ih^Pt = [Ho, pfi. (1) 

Next, part of the theoretical framework of [6] was used to obtain a set 
of linearly independent solutions. Finally, the results were transferred to a 
more general setting. 




Figure 1: Ladder (or cascade) configuration (left), A configuration (center), 
and Vee configuration (right). 



2 Special solutions 

Consider a Hamiltonian of the form 

H = H + acn(Lut,k)[S 4 ] t + xdn(ujt,k)[S 7 ] t . (2) 

where S±, S 2 , ■ ■ ■ , Sg are the generators of SU(3) and equal half the Gell-Mann 
matrices — see the Appendix A — and where 

[Sj] t = e -(it/H)Ho S . e (it/K)E (3) 

are the generators written in the interaction picture. This kind of Hamil- 
tonian is considered in quantum optics when studying three level systems 
driven by laser light, neglecting damping effects — see for instance [U [7]. 

The functions sn , cn , dn are Jacobi's elliptic functions. In the limit 
k = the function sn (ut, k) converges to sin(c<jt), cn (ut, k) converges to 
cos(u)t), and dn (ut, k) converges to 1. In the limit k = 1 the function 
sn (cut, k) converges to tanh(u;t) and cn (cut, k) and dn (tut, k) both converge 
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to sech(ut). In what follows, we drop the arguments (ut,k) of the Jacobi 
functions when this does not lead to ambiguities. 

Let us assume that the parameters of the Hamiltonian satisfy 



4k 2 (fko) 2 = a 2 + k 2 x 2 . 



(4) 



Then three orthonormal solutions ■0o>'0+>'0- °f t ne Schrodinger equation 
ihip = Hip are given by 



ia dn 
k 2 x cn 
Bsn 

(k 2 xT cn ± iaB sn dn 
±k 2 xB sn cn + iaT dn 
Tk 4 x 2 cn 2 T a 2 dn 2 



± e -(it/h)H 



(5) 



(6) 



with 5 = 2k 2 fku and T = Va 2 + A; 4 x 2 . Note that by assumption one has 
T 2 = B 2 + a 2 (I - k 2 ). The functions i?(t) and <p{t) are given by 



V2TVT 2 -B 2 sn 2 = V2T^a 2 {l - k 2 ) + B 2 cn 2 (7) 

1-A; 2 



R(t) : 



k 2 ) + B 2 cn 2 (ujs,k)' 
One verifies the above statements by explicit calculation. 



(8) 



3 Exciting the ground state 

Let us consider a wave function ip(t) which at t — satisfies -0(0) 
It can be decomposed into the basis of special solutions as 

m = -^o(o) + -||(^ + (o)+^-(o)). 

The time-dependent solution is then 



,-(it/h)H 



a 

2^2 



k 2 x 



a dn 

-ik 2 x cn I + V2 
-iB sn / R W 



COs(0(t)) 



+iV2 



k 2 x 



—iaB sn dn 
r> sin(0(t)) I — k 2 xB sn cn 
' 1 - - \ /c 4 a; 2 cn 2 + a 2 dn 2 



(1,0, 0) T 
(9) 




(10) 



Clearly, all three independent solutions are needed to obtain the time evolu- 
tion for the given initial condition. It is also clear that the phase factor e T "^ 
which appears in ([6]) when k ^ 1, although not so relevant for the special 
solutions, becomes highly relevant in the above quantum superposition. 




Figure 2: Occupational probability of level (1) as a function of time for 
k = 0.25, h = lu = 1, a = 0.3 and x = 1.6. 

As expected, the time- dependent interaction populates the two other 
states. See the Figures [2] and [3J Note that level (1) does not go below 
half occupation. 



4 Discussion 

We obtained solitary solutions for a three-level system with periodic time- 
dependent external fields. Two aspects are novel. The external fields are non- 
linear in the sense that Jacobi's elliptic functions are used as deformations 
of the usual harmonic functions. In addition, a full set of special solutions is 
obtained so that arbitrary initial conditions can be considered. 

The additional phase factor exp(=p0(t)) appearing in the solutions (0 [6]) 
was first considered in [6]. It is not very relevant for the special solutions 
themselves, but has effect on their superpositions. The function sin(0(t)) is 
plotted in the Figure HI Its frequency is slightly lower than the frequency 
lu/2tt of the driving field. As a consequence, a low frequency beat appears 
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Figure 3: Occupational probability of level (3) as a function of time for the 
same parameter values as in the previous Figure. 

in the case of a superposition of the special solutions. This is dominantly 
visible in the Figures [2j [3j For the chosen set of parameters the beat period 
is about 5 times the frequency of the external field. Note also the frequency 
doubling by which the exchange of population occurs between levels (2) and 



Figure 4: The function sin(0(t)) for the same parameter values as in the 
previous Figures. 





1.0 — 



5 



Acknowledgments 



We acknowledge discussions with Dr. Maciej Kuna about solving the non- 
linear von Neumann equation in the SU(3)-case. 



Appendix A 

The following expressions are used for the generators of SU(3). 
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& = i ,5, 



S 7 = i -i , S 8 = — ^ 1 . (11 




Appendix B 

The solutions (0 [6]) were obtained starting from a known solution of the 
non-linear von Neumann equation 

3 

ih Pt = -j\{Ho,Pt},Pt], (12) 

where Hq is given by 

/ -fi \ 
F = - y. . (13) 
d \ OA/ 

The three different configurations, vee, ladder, and Lambda, are obtained by 
taking A < —jj, < 0, |A| < respectively A > /x > 0. 

Let H{t) be defined by H(t) = §{# ,p t }. Then p(t) is a solution of the 
linear von Neumann equation with time-dependent Hamiltonian H(t). 

A known solution of the non-linear equation (|T2|) is of the form [51 [8] 

p(t) = -1 + A cn (out, k)[S 4 ] t + Bsn(cut,k)[Si\ t + Cdn(cut,k)[S 7 } t 
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e 



, | ^Bsn |Acn , 

-(it/K)H I l Bsn I _iC dn e (i*/»)Ho j ( 14 ) 



The coefficients A, B, and C, are real. They must satisfy the set of conditions 

frwB = fiAC 
2hiuk 2 C = (X + p)AB 

-2huA = (X-fJ)BC. (15) 

This set of equations can be solved in a straightforward manner when < 
|A| < p. 

Next, a unitary matrix V(t), satisfying 

Pt = v(t) Po v(ty. (i6) 

is calculated. The fastest way to find V(t) is by first diagonalizing p t . Note 
that the eigenvalues of p t do not depend on time. They are given by 

- and -T-T (17) 
3 3 2 y J 



with T = V A 2 + C 2 . The result is 

V(t) = e- {it/K)Ho G{t)G(^ 



with 



G(t) ' 



_ V2TR(t) 

iV2CR(t) dn - AT cn - iBC sn dn AT cn - iBC sn dn 

-y/2AR(t) cn ,45 sn cn + iCT dn AB sn cn - iCTdn | (19) 
y/2BR(t)sa R{t) 2 R{t) 2 



and 



R(t) = VA 2 cn 2 + C 2 dn 2 = VT 2 -B 2 sn 2 . (20) 



However, V(t) does not necessarily describe the unitary time evolution 
U(t). But the latter can be related to V(t) by the method of [B]. The knowl- 
edge of £/(i) implies the time evolution of wavefunctions ip(t) for arbitrary 
initial conditions ipit) = U(t)i/;(Q). It turns out that the special solutions (0 
[6]) are the columns of the matrix G(t), taken in the interaction picture. Two 
of the three solutions are multiplied with the time-dependent phase factor 
exp(=p0(t)). Finally, note that the conditions (1T51) are needed during the 
above derivation but are not required for (jSl E]) to hold. They rather are 
replaced by the single condition (j4j). 
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